arXiv:1505.02882vl [gr-qc] 12 May 2015 


Mon. Not. R. Astron. Soc. 000, fll fT^ (2002) Printed 13 May 2015 (MN style file v2.2) 

Geodesic properties in terms of multipole moments in 
scalar-tensor theories of gravity 

George Pappas^* and Thomas P. Sotiriou^’^’^f 

^ School of Mathematical Sciences, The University of Nottingham, University Park, Nottingham NG7 2RD, UK 
^School of Physics and Astronomy, The University of Nottingham, University Park, Nottingham NG7 2RD, UK 
^Perimeter Institute for Theoretical Physics, Waterloo, Ontatio, N2L 2Y5, Ganada 


ABSTRACT 

The formalism for describing a metric and the corresponding scalar in terms of multi¬ 
pole moments has recently been developed for scalar-tensor theories. We take advan¬ 
tage of this formalism in order to obtain expressions for the observables that charac¬ 
terise geodesics in terms of the moments. These expressions provide some insight into 
how the structure of a scalarized compact object affects observables. They can also 
be used to understand how deviations from general relativity are imprinted on the 
observables. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are astrophysical sys¬ 
tems that host black holes (BHs) or neutron stars (NSs) 
with a regular star as a companion. The presence of a com¬ 
pact objects in the binary makes it a natural laboratory for 
some of the most extreme aspects of physics. In particular, 
the systems that harbor NSs involve both strong gravity ef¬ 
fects and effects associated to the properties of matter in 
supra-nuclear densities. Therefore, studying these systems 
can broaden both our understanding of gravity, by testing 
the predictions of the established theory of General Relativ¬ 
ity (GR), as well as our understanding of the equation of 
state (EoS) of matter for densities as high as the densities 
that can be found at the centre of NSs. 

There are several astrophysical observables that can 
be associated to the properties of the spacetime and more 
specifically to the properties of geodesics around the com¬ 
pact object that is part of such a system. For example, quasi- 
periodic oscillations (QPOs) of the X-ray flux observ e d from 
the accretion disc in LMXBs (for reviews see iLam 3 ll2003ll : 
Ivan der KlisI (I2OO6I II is one such phenomenon. The proposed 
mechanism for expla i ning i t is the relativistic precession of 
geodesics fsee IStellal (I2OOII IL In this context, one assumes 
that the observed QPOs come from the frequencies asso¬ 
ciated to circular geodesic orbits and their perturbations. 
The main frequencies related to these orbits are: the orbital 
frequency of the circular motion fl; the radial oscillation of 
a slightly elliptic orbit that has frequency Kp and the cor- 
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responding precession frequency flp = Q — Kp (periastron 
precession); and finally the vertical oscillation of a slightly 
off-equatorial orbit which has frequency and the corre¬ 
sponding precession frequency = fl — (nodal preces¬ 
sion). As is customary in the relevant literature, the term 
“frequency” is used loosely to refer to what is actually an an¬ 
gular velocity that is equal to 2Tn/, where ly is the frequency 
expressed as cycles in the unit of time. 

The radiated spectrum of an accretion disc can be also 
associated to the properties of the geodesics around a com¬ 
pact object. It is associated to the energy per unit mass E/yL 
of a particle that has a circular equatorial orbit and how that 
is distributed from the outer parts of the disc down to the lo¬ 
cation of the innermost stable circular orbit (ISCO), which is 
essentially the position where the accretion disc terminates. 
Under the assumption of a thin, radiatively-efficient accre¬ 
tion disc, a fluid element of unit mass in the disc radiates as 
it moves inwards in a way that depends on the distribution 
of energy on the circular geodesics. This information is im¬ 
printed on the temperature distribution of the disc and its 
spectrum. This type of properties of an accretion disc and 
their connection to the background geometry have been used 
in the p ast to determ i ne the rotation of black holes (see for 
example IShafee et al.l l|2006ll and references therein), where 
a technique called continuum-fitting method is applied in 
order to determine the spin of the black hole from the lumi¬ 
nosity of the disc. 

One would like to extract information about the struc¬ 
ture of a compact object from observables by studying the 
properties of the geodesics around it. One way to do so is to 
explore the relation between these observables and the mul¬ 
tipole moments of the spacetime as they were defined by 
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iGerochl lll970al f9): iHanseJ lll974) : iFodor et alj (Il989l ) . This 
has been the topic of extensive study in GR. iRvan (Il995li 
studied the relation between the orbital and precession fre¬ 
quencies and the multipole moments of the spacetime. He 
also studied the relation between the multipole moments 
and the energy per unit mass, E/fj., distribution of circular 
geodesics. In particular, he expressed the energy change per 
logarithmic orbital frequency change. 


AE _ 1 dE 

fM fM din D, 


( 1 ) 


as well as the ratios of the precession frequencies over the 
orbital frequency, Qp/Q and as functions of the or¬ 

bital frequency H in the form of series expansions where 
the coefficients of the expansions depend on the multipole 
moments. 

These relations were intended to be used in gravita¬ 
tional waves detection in order to extract the multipole 
moments of the spacetime around supermassive BHs, but 
they have proven to be useful in other ways as well. They 
we re used to calculate the mom ents of numerical spacetimes 
by Laarakkers fc Poisson! (Il999l l and IPapoas fc ApostolatosI 
(l2012tl and it seems that the method has become a common 
practice after that (with the most recent examples the ap¬ 
plication to anisotropic stars by YagL&TYungg (l2015l i and to 
tidally deformed objects bv IPani et al.l ( 2015 1'). They were 
also proposed as a possible way of measurin g the fir s t thre e 
multipole moments of a NS from QPOs by |Pappli3 ll2012ll . 
They h ave even been u sed t o study floating orbits around 
BH s by Kapadiai.gLay 1 2013ll . 

IShibata Sz Sasakil Jl998l i have studied how the ISCO 
radius, i?isco, and the orbital frequency, flisco, at the lo¬ 
cation of the ISCO, depend on the multipole moments of the 
spacetime, with the aim to impose constrains on the prop¬ 
erties of the central compact object from measurements of 
these quantities. 

There is no doubt, therefore, that this sort of rela¬ 
tions between orbital properties and multipole moments 
have been a very useful tool in studying compact objects 
in GR. As a next step, one could attempt to derive similar 
relations for alternative theories of gravity. The motivation 
for this is actually 2-fold. One can obviously use such rela¬ 
tion in the same applications that they have already been 
used in GR. In addition to this though, one could also ex¬ 
ploit them in order to understand the degeneracies in the 
effect that deviations from GR and uncertainties in the de¬ 
scription of the matter content of the star itself can have on 
actual observables. 

Recently, the multipole moments for scalar-tensor the¬ 
ories of gravity with a mas s less s calar field have been de¬ 
fined bv IPappas fc Sotirioul (l2015l ). This opens the way for 
relating geodesic properties that are associated to astrophys- 
ical observables, such as frequencies, accretion spectra, and 
ISCO radii, to the multipole moments in scalar-tensor the¬ 
ory. These relations and the way that they compare to the 
corresponding GR relations, can be a very useful tool for 
testing GR against observations and constraining the pa¬ 
rameters of scalar-tensor theory. 

In what follows, we first give in Sec.[2]a brief description 
of scalar-tensor theory with a massless scalar field, a descrip¬ 
tion of stationary axisymmetric spacetimes in this theory, 
and the properties of circular geodesics in these spacetimes. 


In Sec. Owe calculate expressions for AE/fi, the various fre¬ 
quencies, and the ISCO radius in terms of the moments for 
the spherically symmetric solution in scalar-tensor theory. 
This serves as a first demonstration of the methods that we 
then use in Sec. |4] in order to calculate the same expres¬ 
sions for axisymmetric solutions. In Sec. O we discuss how 
the expressions could be used to distinguish between GR 
and scalar-tensor theory. Finally, in Sec. [6] we present our 
conclusions. In the Appendix we give the full expressions 
for the ISCO radii and the orbital frequencies. 


2 CIRCULAR EQUATORIAL GEODESICS IN 
SCALAR-TENSOR GRAVITY 


One of the most thoroughly studied alternatives 


CaDOzziello & Faraonil 

201 ll: iDamour & Esoosito-Farese 

1992 

Jordanlll949l: Fierz 

19561: 1 Jordan 19591: Brans & Dicke 


1961 


Dickd 1962ll . which can be described by the following 


action 


S = J d'^xV^ , 

( 2 ) 

where g is the determinant and R is the Ricci scalar of the 
metric g^^, denotes the corresponding covariant deriva¬ 
tive, Sm is the matter Lagrangian, and ip collectively denotes 
the matter fields. It is assumed that matter helds couple 
minimally to gp,^. 

In action m the scalar field <1? is nonminimally coupled 
to gravity and has a noncanonical kinetic term. This rep¬ 
resentation of the theory is called the Jordan frame. The 
conformal transformation = 16nG $ g^v, together with 
the scalar held redehnition 


d(p 


2uj{^) -b 3 


dln$, 


will bring action m into the following form 


( 3 ) 


where the matter helds still couple minimally to g^^- This 
implies that (p is now coupled to the matter helds, and it is 
this coupling that encodes any deviation from standard GR 
with a minimally coupled scalar held. This frame is called 
the Einstein frame. The advantage of the Einstein frame is 
that the held equations outside the matter sources take the 
form. 


Rab = 2da(pdb4>, 

r'’^aVb<P = 0 , 


( 5 ) 

(6) 


which are essentially GR with a minimally coupled scalar 
held. 

A spacetime is stationary and axisymmetric if it ad¬ 
mits a timelike Killing vector held G and a spacelike Killing 
vector held 77 “ that has closed integral curves. The actions 
of these symmetries should also commute, i.e., 77 “— 
^“Va77^ = 0. Furthermore, the condition for the 2-planes 
that are orth ogonal to the two Ki ll ing vectors to be inte- 
grable is (see lStephani et ^ ll2003ll : IWaldl (ligSdh l: 

eabcdrj ^ V ^ = 0 = Eabcd^ rj V 7] . (7) 
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This condition can be also written in terms of the Ricci 
tensor in the form: 

Rd[aib'nc] = 0 = rf Rd[aGVc\- ( 8 ) 

In GR this condition is satisfied in vacuum. As a con¬ 
sequence, the line element of stationary and axisymmet- 
ric vacuum spacet imes can take the Weyl-Papapetrou form 
dPapapetroulllOS^ . 

ds^ = — f {dt — wd(p)^ + f~^ (dp^ + dz^) + p^d<p^] , (9) 

without loss of generality, where the metric functions depend 
only on the coordinates (p, z). 

In the case of scalar-tensor theories, one can use the 
vacuum field equations in the Einstein frame, eq. m , in or¬ 
der to show that condition m is satisfied. By virtue of the 
assumption that the scalar field obeys the symmetries of 
the metric, i.e., ^°'Va<j> = 0 and r]°‘Vafj) = 0 , one has that 
GRab = v'^Rab = 0. This implies that the conditions for 
integrability are satisfied and the line element for a station¬ 
ary and axisymmetric spacetime in scalar tensor theory can 
be written in the Weyl-Papapetrou form, eq. m, without 
any loss of generali t y. Int erestingly enough, as it is shown in 
IPappas fc Sotiri^ (1201511 one can write the field equations 
for the Weyl-Papapetrou metric in the Einstein frame in the 
same form as they are in GR. The field equations can, there- 
fore, be writ ten in the form of the equation introduced by 
lErnstl (Il968ll 

(7^^:= Vf ■ Vf, ( 10 ) 

where £ = f + iu}, oj is the scalar twist of the timelike Killing 
yector (not to be confused with the tj(<I>) in the action 
of the scalar-tensor theory), and V is the gradient in flat 
cylindrical coordinates. The Ernst equation is accompanied 
by two equations for the function 7 (which we will suppress 
here for breyity) and one more equation for the scalar field, 

VV = 0 , ( 11 ) 

and this closes the set of equations that characterise a so¬ 
lution of the scalar-tensor theory in the Einstein frame, the 
solution being a spacetime metric Qab and a scalar field 4>. 

If we are to study the properties of geodesics that mat¬ 
ter follows, we need to go from the metric in the Einstein 
frame to the metric in the Jordan frame, which is the one 
with which matter couples minimally. The metric in the 
Jordan frame can be expressed as the metric in the Ein¬ 
stein frame times a conformal factor that depends on the 
scalar field on the Einstein frame, gab ~ (167rG$)“^Pa6 = 
Ad{(j))ga.b, where the Einstein frame metric is expressed in 
the Weyl-Papapetrou form, eq. ®. The conformal factor 
A^((()) is a free function of the theory which can be deter¬ 
mined by integrating eq. for a giyen cj($). 

Since the scalar field is assumed to respect the sym¬ 
metries of the metric, the Jordan and the Einstein frame 
metrics will share the same symmetries. Therefore, in both 
frames there exist two killing fields which define the symme¬ 
tries with respect to translations in time and rotations with 
respect to an axis. One can take advantage of these sym¬ 
metries to study geodesics in the spacetime in the Jordan 
frame, as one would in GR. 

Symmetry under time translations is associated to an 


integral of motion, the energy E 


E = -paG = -pt=m : (12) 

where t is the coordinate time and r is proper time. Sym¬ 
metry under rotations is again associated to an integral of 
motion, the angular momentum L 

L = PaV'^ = Pv = ^ (^9t<p^ + 9v<P^'^ ( 13 ) 

In addition to the two previous equations, from the four- 
momentum of a particle, which is defined as = mu°' = 
we have the equation 


- 1 = 9tt 


dt\^ f dt\ f dip 


. dtp 

+ 1 ^ 




dz 


(14) 


If we further define the angular velocity, ft = ^, for the 
circular and equatorial orbits, eq. (US defines the redshift 
factor between coordinate and proper time, 


) — 9tt , ( 13 ) 


and the energy and the angular momentum for the circular 
orbits take the form, 




-gtt - gt<p^ 


g -\J gtt ‘^gtip^ g<p<p^^ 


9 


9tip + gip^p£i‘ 


\J gtt ‘2gtip£i g^p^p£t^ 


(16) 

(17) 


where we have introduced the energy and angular momen- 

, .2 

turn per unit mass. From the conditions, -jf = 0, = 0 

and ^ = 0 for circular equatorial orbits, and the equations 
of motion obtained assuming the Lagrangian, C — ^gabi°‘T’, 
the angular velocity can be calculated to be. 


Q. = 


gtip,p + G {gtip,pT 9tt,p9w,p 
9'p‘p,p 


(18) 


This is the orbital frequency of a particle in a circular orbit 
on the equatorial plane. 

Equation m can take a more general form in terms of 
the constants of motion, 


dp 


dz\ _ ^ E g^pip 2ELgtip L gtt 
dr) {gtvY - 9tt9vv 


= Vaff. 


(19) 


With equation GSI we can study the general properties of 
the motion of a particle from the properties of the effective 
potential. Furthermore, we can study perturbations around 
circular equatorial orbits. If we assume small deviations of 
the form, p = pc + 5p and 2 ; = Sz, then we obtain the 
perturbed form of m, 


(d[5p)V 


d{Sz) 

dr 


1 d^Vcff 

2 

IdWcff 


+ 2 


dz^ 


iSpf 

(Szf. 


( 20 ) 
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This equation describes two harmonic oscillators with fre¬ 
quencies, 


- 2 _ 

^ 2 C»p2 

-2_ gVe// 
" 2 dz^ 


( 21 ) 

( 22 ) 


The differences of these frequencies (corrected with the red- 
shift factor m, i-e., Ka = [dr/dt)Ka) from the orbital fre¬ 
quency, Q.a = O. — Ka, define the precession frequencies, 
where the oscillation frequencies Ka are given in terms of 
the metric functions as. 


-2(gu + gt^mgt^ + g^^^) (23) 

\2 / 9tt 




\A'^{4>)p^ 


where the index a takes either the value p or z to express the 
frequency of the radial or the vertical perturbation respec¬ 
tively and the expressions are evaluated on the equatorial 
plane z = 0. The position where k^ becomes zero is the 
location of the ISCO. These equations will be used to calcu¬ 
late -Risco, the various frequencies flp, flz and fl, and the 
quantity AE = — in scalar-tensor theory. 

There is one final issue that needs to be discussed. Prop¬ 
erties associated with geodesics of a metric will depend only 
on the metric in question and not on the corresponding 
scalar field. Here we are considering geodesics of the Jordan 
frame metric, as these are the ones followed by test particles. 
Their properties will not depend directly on >I> (obviously $ 
does implicitly affect the geodesics by acting as a source in 
the modified Einstein equations). One can be tricked into 
thinking that the scalar moments will therefore not appear 
in the expression for the various observables, as the met¬ 
ric is fully determined by its own moments. Ho wever, the 
multipole moments in IPappas fc Sotiri^ (l2015l l are actu¬ 
ally defined in the Einstein frame and ( 167 rG$)“^ = 
is the conformal factor that relates the Jordan and the Ein¬ 
stein metrics. So, the observables will explicitly depend on 
both the scalar and metric multipole moments. 

The hidden presence of the free function A{(j>) of the 
conformal factor also highlights the need to select a specific 
function A{(j>) in order to pin down the theory in consider¬ 
ation [recall that A{(fi) is related to i^(<?!>) through eq. ([3)l]. 
We would like to avoid doing so in order to keep our results 
as general as possible, so we will adopt the following expan¬ 
sion of the conformal factor around the asymptotic value of 
scalar field, t^oo = <j>o, 


A{(f>) = A{(j>o) -I- A'{(j>o){(l> - M + ^”('(>o)(<)> - 
+ A® (00) ((() - 0o)^ + • ■ • 

Here, we can define in the usual way (see 


for example 


_ _ Da mour fc Esposito-Faresd lll992ll : 

IPamour fc Esposito-Earesd ( 1996fl l the parameters 
ao = dlnA/d 0 |oo, do = da/d 0 |oo, 7 o = d 0 /d 0 |oo, 

and so on. The derivatives of A{(fi) at infinity can then 
be expressed in terms of the parameters a, 7 , etc. (for 
example, A”{(po) = A((po){d + od))- 

With this parameterisation of the conformal factor any 
observable associated to geodesics in the Jordan frame can 
be expressed in terms of the moments and the coefficients 
of the expansion, i.e., A{(j>o), q,/3 ,7, etc. Note that the 
asymptotic value of the Jordan frame scalar <l?oo is related 
to the asymptotic value of the Einstein frame scalar 0o 
via (IOttG'I’ oo)”^ = H^(0o)- In what follows, we will set 
A{(j)o) = 1 for simplicity. This amounts to a Weyl rescal¬ 
ing of the Jordan frame metric and it should be taken into 
consideration when one interprets the final expressions for 
the observables. In particular, according to their definitions 
n, Qp, and flz are insensitive to such a rescaling, whereas 
AE/p, picks up an A{(f>o)~^ factor if p is held fixed, p actu¬ 
ally corresponds to some particle mass, p^ = —PaP“, so in 
reality this mass would be also rescaled by a factor A{(/>o)~^ 
and the ratio AE/p would be insensitive to constant Weyl 
rescalings. Caution should also be taken when calculating 
Q, d, 7 , etc. for a given theory, as their values are affected 
by such rescalings. 

Finally, we should stressed that the multipole moments 
for the scalar field are defined up to a constant shift of the 
Einstein scalar, i.e. they actually characterise the quantity 
d — <j>o- When calculating the moments one can set 0o = 0 
without loss of generality and this is the approach we will 
follow in the rest of the paper. As is obvious from the above 
though, the actual value of do does affect the observables 
through A{do) and the higher order coefficients of the ex¬ 
pansion for A{d) and this will be taken consistently into 
account. 


3 OBSERVABLES AND MOMENTS IN 

SPHERICALLY SYMMETRIC SOLUTIONS 


As a warm-up calculation, we will start with the static, 
spherically symmetric solution in scalar-tensor theory. The 
metric, expressed in the Einstein frame, is 


ds^ 



dP + 



-k 



P{d9^ + sin^ 9difi^). 


(25) 


In these coordinates, the corresponding scalar field is 


<^=^log (1 



(26) 


and rrP + w\ — P, i.e. this is actually a 2 -parameter family 
of solutions. This is the static, spherically s ymme t ric so ¬ 
lution in scalar-tensor theory tha t is given in ljus 3 lll959ll : 
iDamour fc Esposito-Fare^ (ll992h P1 


^ In the notation of lPamour fc Esposito-Faresd lll992h 21 corre- 
sponds to a, 2m, corresponds b, and the scalar charge wa corre¬ 
sponds to d. The bond m?-\-w'^ = corresponds to a^ — b^ = Ad?. 
The comparison also reveals that m is the Einstein frame mass. 
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The metric written in Weyl-Papapetrou coordinates 
takes the form 


ds = — 


where 


R-+R+- 21 
R- + R+ + 21 


1 I z^+p-l^ 

2 2R+R- 


dT + 


R— + — 21 

i?_ + i?+ + 21 

{dp^ + dz^) + dG^ 


R+ = ^{1-zY+P^., 


R- = G(l + zf+p^. 

The scalar profile in the same coordinates is 


WA , 


Gd + zY + p^-i-z 
Gd-zY + p^ + i-z 


(27) 

(28) 

(29) 

(30) 


where wa turns ont to be the scalar charge. 

In order to calculate the orbital parameters that one 
would measure for test particles moving along geodesics in 
this spacetime, one needs to use the metric in the Jordan 
frame, which is conformally related to the metric in the 
Einstein frame through the relation gab = A^{(j))gab- In this 
expression then, one can substitute A{(j}) with its Taylor ex¬ 
pansion in powers of — 4>o), which can then be replaced 
by eq. (l3?il) . 

Following iRvanI (Il995ll we would like to express the var¬ 
ious observables as series expansions in the orbital frequency 
fl. The coefficients of such an expansion will end up depend¬ 
ing on the multipole moments and one could in principle 
then determine the latter order-by-order from observations. 
The orbital frequency of circular orbits on the equatorial 
plane (z = 0) can be shown to have the following expansion 
in inverse powers of p, 

M + WACtO^ ^ I . . -l/2\ /Qi\ 

-g- 1 (^1 -I- series m p ' j . (31) 

We can see from the expansion that we can redefine the mass 
as 



M = M + waolq ■ (32) 

It is preferable to work in terms of the dimensionless quan¬ 
tities U = (MII)^^® and x = {MjpY^'^■ U can be expressed 
as a series expansion in x as 


2qoM — PoWA 

[2w'’AF(37o - 2ao/3o) 

+ 2w'^M^ (llal+gPo-S) 

-2AaowM^ + 9M* - 5,5ow'‘] + O {x^ (33) 

This series expansion can be inverted and x can be expressed 
as an expansion in U. Observables can then be expressed as 
power series in U. Note that U is preferable to Mfl because 
the expansion will contain integer powers of U, as opposed 
to fractional powers of Mtl. 

With the strategy laid out, we can proceed to deriving 
the actual expressions for the various observables, starting 



U=x+- 

6 


WA 


with AE/p. The energy per unit mass of circular orbits is 
given by eq. Gi). Once this is expressed as a series expansion 
in U, the energy per logarithmic frequency interval is given 
as 


AE _ U_GE 
~ ^ ~JdU’ 


(34) 


which as a series has the following expansion. 


17= 

r ‘^Pow'a 

8aoWA 1 

rr4 , 


9M2 

9M 2 

u + 


8Pow\ 


— 8MwA(4ao/3o + 7o) + iw^M^ (2 — 38ao + 3/9o) 


-h 22‘iaoWAM^ - 81M'‘ 




2AM* 


+ O (I/®) . (35) 


If one sets the scalar field charge to zero, one recovers the 
Schwarzschild expansion in GR. It is already evident in this 
expression that one cannot do away with the effects of the 
corrections coming from the scalar field by a redefinition of 
the mass. 

One can proceed to calculate the expansions for the two 
epicyclic frequencies, tip and tlz- The periastron precession 
will give the expansion. 


np 

n 


3- 


WA (saoM -\- Powa'^ 


2M2 




4m(2w\M {I7al + 3/3o - 3) - GOagWAM^ 
27M® -h 37owi) - 13PWa ^ 


2AM* 


APow\M{2aoPo + II 70 ) - (20qo;5o + 6ao7o 


Pq -i- APq -h 2(5o) + SwaM^ ^ — 112rrQ -t- 2ao(/3o + 8) 

-h 370 ) -h 12w\M^ (l96ao -h 7Po - 12) - 2112aoWAM® 


rjo 

+ 6A8M’'- 35P^w\ -S^ + o{uG 
48M6 ^ ^ 


(36) 


Again, this expression reduces to the Schwarzschild one if 
we set the scalar charge equal to zero. 

As is the case in GR, the spherically symmetric solution 
in scalar-tensor leads to zero nodal precession, since Kz = tl 
and, therefore, tlz/tl = 00 

What remains to be calculated, is the location of the 
ISGO. One way of identifying the location of the ISCO is 
as the radius at which the radial oscillation frequency of a 
slightly eccentric orbit becomes zero. Therefore one needs to 
calculate the roots of the square of the frequency Kp. Even 
in the simple case of spherical symmetry, the evaluation of 
the roots is not easy and there is additionally the problem 
of having an unspecified function, which is the conformal 
factor for the transformation to the Jordan frame, A{(f>). 


^ In the work bv iDoneva et all ||2014| ~). one can find analytic ex¬ 
pressions for the orbital frequency and the radial oscillation fre¬ 
quency of the spherically symmetric spacetime in scalar-tensor 
theory, for a specific choice of the conformal factor. 
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This problem can be solved using a perturbative approach 
as long as the compact object is not considered to be too 
scalarized. In that case, one could treat the scalar charge 
(actually the scalar charge over the mass) as a perturbative 
quantity of order e, while at the same time assume that the 
sought after root has the GR value plus corrections, i.e., 
Pisco = 2\/6M{l + cie + C 2 £^ + ...)• This expression for 
the radius can be substituted in the function Kp, which can 
then be expanded as a series in e. The result will be a sys¬ 
tem of equations that will have to be solved in terms of the 
coefficients Ci so that the expansion gives zero at all orders. 
T his calculation is analogo us to the calculation presented 
bv IShibata fc Sasakil (Il998h . The resulting expression up to 
fourth order for the ISCO is. 


4 OBSERVABLES AND MOMENTS IN 
AXISYMMETRIC SOLUTIONS 

We will now derive the general expressions that relate the 
orbital frequency to the precession frequencies and the en¬ 
ergy change per logarithmic frequency interval for a station¬ 
ary and axisymmetric spacetime in scalar-tensor theory of 
gravity. We will also present expressions for the location of 
the ISCO and the orbital frequency at the ISCO which will 
hold under certain reasonable assumptions regarding the be¬ 
haviour of the higher order moments of both the spacetime 
and the scalar field. 


4.1 Ernst potential, scalar field and moments 


Pisco =2V6m(i - + (0.0371084ag - 0.092897/90 

2 

-0.072918) ^ -t [o.009646ao + ao(0.197845/3o 

3 

- 0.0201845) -f 0 . 02525570 ] ^ -t [0.000736 qo 

-f ag(-0.185728/3o - 0.0000175) - 0.0290932ao7o 
-f (-0.0225927/30 - 0.001873)/3o - 0.0027775o 

4 

- 0.0046151] ^ + ■ ■ ■ ) (37) 


which in terms of circumferential radius can be written as, 

0.452733aoWA , nonncoc 2 
-Risco —6A/ --|- (^0.0399826ctQ 

2 

-0.053767/9o - 0.059701) ^ + (0.0115299ao 

3 

-tao(0.159059/3o - 0.017441) + O.OI 88 I 570 ) 

+ [ao(0.0001688 - 0.143397/3o) + 0.001467ao 
- 0.0232982ao7o - (0.0154855/3o -t 0.000099)/3o 

4 

-O.OO 21512 J 0 - 0.003793] ^ + ■ ■ • ) (38) 

Finally, one last useful quantity at the ISCO that we could 
evaluate in terms of the parameters of the spacetime is the 
orbital frequency, 

IIisco =—^ (1 -t + (1.21083 X 10-®Qo 

gVgm V M '' 

2 

-f0.0354517/3o + 0.0777343) ^ + (ao [0.089554 

3 

-0.117187/3o] +0.007811qo - 0.0225997o) ^ 

-f (-0.0424804qo + ao(0.009887/3o +0.023437) 

+ 0.009229ao7o + (0.0178996/3o + 0.006044)/3o 


4 

Wa 


+0.0028125o + 0.00983) ^ + 




(39) 


One could ask whether the assumption of a weakly 
scalarized compact object is a reasonable one. The current 
constrains on the values of the parameters qq , /3o for scalar- 
tensor theories are such that do not allow for very large 
scalarization for non-rotating compact objects. Therefore, 
treating the scalar field as a perturbative parameter in the 
non -rota t ing ca s e is a reason a ble a ss umption [for exampl e 
see IWilll ll200(^: iFreire et al.l ll2012li: iDoneva et al.l ll2013li : 
IPani fc Bertil (120141: Berti et al.l (2015h]. 


As we have briefly discussed in Section [2] in the Einstein 
frame the stationary and axisymmetric solutions of scalar- 
tensor theory for a massless scalar field, can be described in 
terms of an Ernst pot ential and a scalar f i eld. In particular, 
as it is described by IPappas fc Sotirioul ll2015li , the Ernst 
potential and the scalar field of the solution are fully de¬ 
termined by their values along the axis of symmetry, which 
can in turn be given in terms of the multipole moments 
of both the scalar field and the metric. We remind that 
in the calculation of the moments the auxiliary potential 
^ = (1 — S)/{1 + £) is also used. 

We define the coordinates at infinity p = p/(p^ + P), 
z = 2 /(p^+ 2 ^) where P = -\-P. The potential f = (1/r)^ 

can be expressed as a series expansion around infinity of the 
form, 

00 

1= ^ aijp'P, (40) 

i,j=0 

where the coefficients aij can be expressed with respect 
to the coefficients aoj = rtij oi the expansion of ^ along 
the axis of symmetry, ^(p = 0) = ■ The coeffi¬ 

cients rrij are related to the multipole moments of the space- 
time and can be found b y solving the expressions given by 
IPappas fc Sotiri^ (l2015l i. In a similar way, one can express 
(j) = {l/f)(f> as a series expansion at infinity of the form, 

00 

^ = ^bijp’-p, (41) 

i,j=0 

and one can see that the coefficients bij can be calculated in 
terms of the coefficients boj = Wj of the expansion of (j> along 
the symmetry axis, (j){p = 0) = '^JLg'WjZ^. The Wj coeffi¬ 
cients are related to the scalar field mom ents and can be also 
calcul ated from the expressions given bv IPappas fc Sotirioul 

iioTi). 

Once we have the Ernst potential and the scalar field, 


^ r - i{p,z) ^ Kp,z) 
r + I(P, z) ’ r 


(42) 


where r = pp + P, we can proceed to calculate the metric 
functions. 

The function /(p, z) of the metric ((9]) is the real part 
of the Ernst potential / = K(f’), while the imaginary part 
of the Ernst potential is cj = Q'(f)- This uj{p,z) function is 
related to the metric function ui(p, 2 ), in ((9l), by the identity 
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Geodesics and moments in ST gravity 7 


/ = —p X Vw, (43) 

where V is the gradient in the cylindrical flat coordinates 
[p,z,tp) and n is a unit vector in the azimuthal direction. 
This equation can be integrated to give the metric function 
ui(p, z) in terms of the moments (for example see the dis¬ 
cussion by iRvarJ lll995ll i. The only metric function that re¬ 
mains to be determined is the function 'y (p,z) and it can be 
evalua ted by integrating eqs. (47)-(48) in IPappafi fc Sotirioul 
(l2015h . in the same way as one w ould do in GR (as before, 
see the discussion by iRva 7l (ll995h l. 

This calculation will give us the Einstein frame met¬ 
ric Qab in terms of the moments. The Jordan frame metric, 
which is the metric to which particles and photons couple, 
will then be given as Qab = A^{(f))gab- This metric is now 
expressed in terms of the spacetime multipole moments and 
the moments of the scalar field, and is the metric that we 
will use to calculate all the quantities that are related to the 
geodesics. 

We should note here, that we assume that particles and 
photons follow the geodesics of the Jordan frame metric and 
that there are no interactions of the particles and photons 
with the scalar field. If we were to assume such interactions, 
then the orbits would deviate from being geodesic. 


4.2 Energy and frequencies of equatorial orbits 

In what follows, we will adopt the following convention for 
the multipole moments of the spacetime and the moments 
of the scalar held in c onnec tion to the moments defined 
by IPappas fc Sotirioul ll2015ll : the mass moments will be 
Ma = K(Pa)! where in particular we will have for the mass 
Mo = ?R.{Pq) = mo = M, the angular momentum moments 
will be Ja = Q{Pi), where the angular momentum will be 
Ji = Q{Pf) = £y(mi), and finally for the scalar field mo¬ 
ments we will have Wa = P^, where the scalar monopole 
will specifically be Wo = Pq = u>o- 

We proceed as we did in the spherical case by evaluating 
the orbital frequency of circular equatorial orbits. In this 
case the expansion of the orbital frequency in inverse powers 
of p will be 

^ (M -Woao^'^^ , . . 

Q, = I -g- I (^1 + senes m p ' j . (44) 

Here we redefine the mass as M = M — Wooto, noting that 
there is a difference in the sign convention for the scalar 
field that results in having essentially Wo = —wa- Again, we 
expand U = as a series in powers of a: = {M/p)^^^ 

and invert the expansion so that we may express a; as a series 
in powers of U. From eqs. m and (1341) we finally get to the 
expansion of the energy change per logarithmic frequency 
interval change in terms of U, 


AE _U‘^ /2/3oWo SaoWo 1\ 4 20Ji[/® 
p ^ 3 V 9M 2 J 9M2 

-l- ^8M [ 3 M 2 + ITo (4 q!o/3o + 7o) ~ 300 W 2 ] 

-b 4Wo (-38ao -b 3/3o + 2) - 224aoWoM^ 

28Ji (lOaoWoM + 9M^ -b 2/3oVFo) 7 / g, 

27M3 '' '' ■ 

(45) 

This expression reduces to the one corresponding to the non¬ 
rotating case of the previous section if we set to zero the 
higher moments and choose the scalar monopole to be Wo = 
—Wa as we have already mentioned. 

In the same way, we can proceed to calculate the expan¬ 
sions related to the two precession frequencies, i.e., t}p/Q 
and flz/fl. These expressions are 

Hp Wo {PoWo - 8 qoM) \ 2 

n ^ " 2M2 j 

-b ^ — V2M [8M2 -b 7 oIRo — 300IR2) 240ctoW^oAf^ 

+ 8WSM^ {I7al + 3po - 3) + 108M^ - l3plW^) 

2 [^Ti (lloioWoM -b 15Af^ -b 5/IoIFo)] -,.5 

3M4 

-b ^ — QM^ [^16Ti -b Wo (24aoAf2 + 20oo/3oIFo 

-b(6Q!o7o + /3o + 4/3o + 2So)Wo — 24a!oIR2 — 24poW2)\ 

- 8M^ [63Af2 + Wo (—112 qq -b 2 olo{Po + 8) -b 870) 

—9a!oIF2] — iPoWoM [33M2 -b Wo (2aoPo + II70) 
-33aoIT2] -b I2IT0 (l96ao -b 7/3o - 12) 

+ 2112qoIToM® -b 648M® - 35/3^ITo® j 

- ^SOJiM [15M2 -b WoicnoPo + 570) — 15aoIF2] 

-b [2J1W0 (5Wo [376qo + 33/3o -b 39] - 243) 

- 810J3] -b 7980aoJiWoM^ + 4320JiM"* 

+ 455PUiWo^)^ + 0{U^) (46) 

and 

Qz _2JiU^ 3C/'‘(M2 - Q0IF2) 2J1W0 {PoWo - aoM) 

'h' “ M2 ^ 2M^ W 

-b (^M|6M2M + 14Ji^ + Wo [{2al - Wo)W2 

- 2Q0M2] } -b 5/3oIFo (M2 - aoW2)) ^ 

+ ^M^SJi |^33M2 -b Wo (370 — 4qo/3o) — 33aoW2j 

- 2M[30ao JiVFoM -b JiWo (5[5qo - 6] Wo -b 27) 

+ 45 J3] } + 40P^ Ji Wo^) + O ([/«) (47) 

Again we can see here that if we take the angular momentum 
and all the higher order multipole moments to zero, then the 
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8 G. Pappas and T.P. Sotiriou 


nodal precession expression, Slz/n, reduces to zero, while 
the periastron precession expression, ^Ip/Q, reduces to the 
expression for the non-rotating case (taking into account the 
sign reversal in the scalar monopole). 


4.3 Location of the last stable circular orbit 


We now turn to the estimation of the location of the ISCO 
as the locus of points where the radial oscillation frequency 
vanishes. The situation in this case is a little more com¬ 
plicated than what we saw in the spherically symmetric 
case. The perturbative treatment that we followed in the 
spherically symmetric case used the scalar charge over the 
mass as an expansion parameter, e, and made the ansatz 
Pisco = 2y/6M{l -|- ci£ -|- C 2 £^ -l- ...) for the location of the 
ISCO. With this ansatz one can solve equation Kp = 0 or¬ 
der by order in e for the various a coefficients. Once one 
abandons the assumption of spherical symmetry, one needs 
to also account for the rotation. The most accurate way to 
proceed would be to consider two different expansion pa¬ 
rameters, one associated to rotation, where an appropriate 
order in £i would be introduced to all higher order moments, 
and one to scalarization, where an appropriate order in £2 
would be introduced to all scalar moments. 

This would obviously complicate things and would not 
allow for the same treatment as in GR or as we had in the 
spherical case. Therefore, we will prefer to stick to using 
only one expansion parameter and attempt to assign an ap¬ 
propriate order in this parameter e to all higher order mo¬ 
ment, both those of the metric and those of the scalar field. 
Clearly, this will limit significantly the applicability of our 
treatment, but it will provide some insight into how the 
scalar moments affect the location of the ISCO. Note that, 
even though a detailed study of how the various moments 
behave for scalarized stars is still pending, the preliminary 
results that a l ready exist in the literature, e.g. the work by 
iDoneva et al.l ll2013l 'l. can help us make an educated guess 
when assigning an order in e to the various moments. 

Starting with the behaviour of the spacetime moments, 
we will treat the rotation as a perturbation and we will in¬ 
troduce an expansion parameter associated to the angular 
momentum of the star. Therefore, we will assume for the 
spacetime moments that the angular momentum Ji will be 
of order £, the quadrupole moment M 2 will be of order 
the spin octupole J 3 will be of order £®, and the mass hexade- 
capole M 4 will be of order £^. This ass umption is justified by 
t he fa c t that, as it has been shown by iLaarakkers fc PoissonI 
lll999ll : IPapDas fc Apostolatoj (120121 'i nYagi et al.l ( 2014l ~l for 
compact objects, the moments higher than the angular mo¬ 
mentum scale as ji", where j — Ji/M^ is the spin parameter. 
Therefore, if we assume that our perturbative parameter is 
the rotation, then the rest of the moments should scale as the 
corresponding powers of j. Clearly, this behaviour is what 
we see in GR and it is not necessary that it will persist in 
scalar-tensor theory. However, in absence of a better guess 
it should serve as a good first approximation. 

The situation is more complicated for the multipole mo¬ 
ments of the scalar field. As we argued previously, we would 
like to keep one expansion parameter instead of introduc¬ 
ing two, but we have essentially two types of effects, the 
scalarization and the deformation due to rotation. With re¬ 
spect to rotation, iDoneva et al.l (l2013h have shown that it 


introduces some ellipticity in the profile of the scalar field 
with respect to the spherical profile of the non-rotating case. 
Therefore, one expects that the higher order moments of the 
scalar field will receive a contribution that scales with the 
angular momentum in the same way as the corresponding 
mass moments of the spacetime, i.e., the contribution to 
the scalar quadrupole W2 will be of order £^ and the con¬ 
tribution to the scalar hexadecapole W4 will be of order 
e*. If we were to keep the scalar charge as an independent 
expansion parameter, £ 1 , then a reasonable guess for the 
behaviour of the scalar moments would be the following: 
the scalar charge would be Wo ~ £ 1 , the scalar quadrupole 
would be W2 ~ £i£^, and the scalar hexadecapole would be 
IT4 - £l£^ 

In order to maintain one expansion parameter one could 
choose to look at the regime where £ ~ £ 1 . This is the regime 
where the influence of the scalar charge is comparable to 
that of the rotation. In this case, we should assign to the 
scalar monopole an order of Wo ~ £, while the higher order 
moments will scale as, W2 ^ £^ for the scalar quadrupole 
and W4 ~ £® for the scalar hexadecapole. This scaling of the 
scalar moments might seem reasonable, but it is at best just 
an educated guess. One could imagine a situation where at 
a more rapid rotation rate the higher scalar moments scale 
with rotation at a weaker power of e. Therefore an alterna¬ 
tive scaling that one could consider would be to assign to the 
monopole an order of Wo ~ £, but for the higher order mo¬ 
ments to assume that they scale at one order lower than in 
the previous case, having therefore a scalar quadrupole that 
behaves as W2 ~ £^ and a scalar hexadecapole that behaves 
as W4 ~ £^. Of course this is an ad hoc choice, but we will 
use it as an alternative example for illustrative purposes. 

Using the hrst of the two different scalings for the mul¬ 
tipole moments and the usual ansatz for the position of the 
ISCO, i.e., Pisco = 2'\/6AI(l + ci£-|-C 2 £^-|-. ..), the resulting 
expression up to fourth order for the ISCO islj 


Pisco =2Vgm(i + 


5aoWo 


5Ji 




16M 3^M2 

0.092923/3oWo^ 0.072877Wo^ 


+ 


0.037297SalWi 


M 2 

O. 21977 M 2 


M 2 

0.009463agWo® 


+ 


M3 

0.02525970 VCo® 

M3 

O.O68I3IJ3 
M4 

O.lSlSieagJiIUo 

M4 

0.000921a^Wo‘ 

M4 

0.029101qo7oWo‘ 

M4 

0.002777(5oI< 


M 2 

Q.197912aol3oWi 

M3 

0.073633aoIU2 


M3 

0.019801aoIUo 

■ M3 

0.162089aoJiIUo 


M3 
0.289247Ji 


M3 

0.013148aoM2lUo 


M* M* 

0.037135l3oJiWg 0.0678245 JiIUo^ 


M* 

0.185879ai/3oVTo^ 

M4 

0.022607,5glUo^ 
M4 


-b 


M4 

0.000059ailUo^ 

M4 

0.001943/IoHo^ 

Ml 


0.004428IUo 0.094989a§WoIU2 


-b 


Ml 

0.0326777/3oWoIU2 

Ml 


+ 


Ml ^ 
0.005692H/olU2 
Ml 


Ml 


...) (48) 


^ Here we present the expressions up to terms of order M The 
full expressions are given in the Appendix. 
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Geodesics and moments in ST gravity 9 


which in terms of circumferential radius can be written as, 


-Risco =Qm{i + 


0.452733aoVKo 0.544331Ji 0.040133aiWo^ 


M 

0.053788,goWo 0.059668VKo^ 0.110743 qo Ji Wo 

M2 “ M2 “ M3 

O.I79922M2 0.011388agWo® 0.159112ao/3oWo® 


M3 


M3 


M3 


0.017135aoWo O.OI88I970 Wj* 


M3 


M3 


0.054768aoW2 0.054505 J3 0.226244J? 

M3 M4 M4 

0.138986agJiWo 0.011282/3o Ji Wq^ 


M4 

M4 

0.054650Ji Wo 

0.002465aoM2Wo 

M4 

M* 

O.OOlOOOaoWo^ 

0.14352ag/3oWo^ 

M4 

M* 

0.000237aoWo‘ 

0.023304ao7oWo‘ 

M4 

M4 

0.015496/3oWo‘ 

0.000155/3oWo‘ 

M4 

M4 

0.002151(5oWo^ 

0.003643Wo‘ 


M4 


M4 


^ 0.075878agWoW2 0.025303/3oW0W2 


M4 

0.004590Wo W 2 

M4 


M4 




(49) 


Similarly, the orbital frequency at the ISCO is given, up to 
fourth order, by the expression 


flisco — 


6V6M 


1 - ~ \; " + 

/2 A'7'7«2 0'7T;r/2 


M 


M2 


M2 


^ 0.035481/3oWo^ 0.077687Wo2 0.23433M2 


M2 


M2 


M3 


0.509757qoJiWo 0.00784a8Wo^ 0.117236ao/3oWo^ 


M3 


M3 


M3 


+ 


0.089082qoWo® 0.0226047o Wo® 0.065506aoW2 

M3 + aP Xf3 

0.780664J? O.20404IQ0M2W0 0.0425199a^Wo‘ 


M4 


M4 


M4 


^ 0.009988ag/3oWo^ 0.022952ag Wo^ 


M4 


M4 


^ 0.009234ao7oWo^ ^ 0.017918,ggW q^ 0.006122,5oWo^ 


M4 


M4 


M4 


0.002812(5oWo‘ 0.009612Wo 0.037394ag W 0 W 2 

+ M4 ^ M4 M4 

0.033029/3oWoW4 0.006338WoW2 


M4 


M4 


^ 0.114417a§JiWo ^ 0.075325/3oJiWo 


M4 

0.168052Ji Wo 
M4 


Ml 




(50) 


If we now use the second scaling for the multipole mo¬ 
ments, the resulting expression up to fourth order for the 
ISCO is. 




1 + 


5aoWo 5Ji 0.037298agWo® 


16M 3y6M2 


M2 


0.092923/3oWo 

0.072877Wo® 0.21977M2 

M2 

M2 M3 

0.009463agWo® 

0.1979119ao,goWo® 

M3 

M3 

0.019801aoWo® 

0.0252670 Wo® 

M3 

M3 

0.073633aoW2 

0.162089aoJiWo 


-f 


M3 M3 

O.O68I3IJ3 0.289247Ji 0.013148aoM2Wo 


M'l 


M'l 




^ 0.000921agWo^ _ 0.185879ag/3oWo^ 0.000059ag Wq' 


M4 


M4 


M4 


0.029101ao7oWo‘ 0.022607/3o Wq 


M4 M4 

0.001943/3oWo‘ 0.0027775o Wo‘ 


M4 


M4 


0.004428Wo^ 0.094989agWoW2 


M'l 


M-i 


0.032678/3oWoW2 0.005692WoW2 

+ M4 ^ M‘^ 

0.181816a§JiWo 0.037135/3o Ji Wo^ 


M4 

0.067825Ji Wo 
M4 


M^ 


+ 


...) 


(51) 


0.75aoWo 0.748455Ji 0.000218ag Wo^ 


which in terms of circumferential radius can be written as, 

r, 0.452733aoWo 0.544331Ji 

R,,,, =6M(1 +- - - 

0.040133agWo= 0.053788,90 Wq^ 0.059668Wo® 


-f 


M2 M2 M2 

0.110743aoJiWo 0.179922M2 0.011388ag Wq® 


M3 


M3 


M3 


0.159112ao/3oW| 0.017135aoWo® 


-f 


M3 M3 

O.OI 88 I 970 W 0 ® 0.054768aoW2 

M3 + W 

0.054505J3 0.138986agJiWo 0.011282/3o Ji Wq 

M* M* 

0.05465Ji Wo 0.226244Ji® 


M4 


M* 


M* 


0.002465aoM2lW ^ 0.001609aoWo‘ 


M4 


M4 


0.14352ag,9oWo^ 0.000237agW( 


M'l M* 

0.023304ao7oWo‘ 0.015496/3o Wq^ 

M4 M4 

O.OOOlSS/SoWo' 0.002151(5oWo‘ 


M4 


M-i 


0.003643Wo ^ 0.075878agWoW2 


M4 


M4 


^ 0.025303/3oWoW2 0.00459WoW2 


+ ... ) (52) 


M* M* 

Finally, the orbital frequency at the ISCO is given, up to 
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10 G. Pappas and T.P. Sotiriou 


fourth order, by the expression 

o _ 1 0.75aoWo 0.748455Ji 

w— + 

0.000218aoVKo 0.035481,5oH^o 

M2 M2 

0.077688Wo 0.23433M2 0.00784ag VKq® 

M 2 M3 “ M3 

0.117236ao/3oW"o 0.089082ao VKq® 

M3 “ M3 

O.O 226 O 470 VF 0 ® 0.065506aoVK2 

W M3 

0.509757aoJiVKo 0.7806640.076216J 3 

W M4 M4 

0.204041aoM2Wo 0.042520^14^0^ 

M4 M4 

0.009988oo/3oVKo 0.022952o§4Ko^ 

+ M4 M4 

0.009234oo7o14"o^ 0.017918/3o VKo 

+ M4 M4 

0.006122/3ol4/o^ 0.002812(5o4Vo‘ 

+ M4 M4 

0.0096124Vo‘ 0.037394ag44"o4V2 

+ M4 M4 

0.033029/3o14"o14^2 0.0063384Vo4V2 

M4 M4 

0.114417oo^i4Vo^ 0.075325/3oJiVKo 

+ M4 M4 

0.168052 JilVo^ 

It is straightforward to verify that these expression re¬ 
duce to the expressions of the spherical case, once the higher 
order moments of the spacetime and of the scalar held are 
set to zero (modulo variations due to numerical accuracy). 


+ 




(53) 


5 COMPARING GR TO SCALAR-TENSOR 
GRAVITY 


The comparison of the expressions that we have produced 
here with the corresponding expressions in GR, lead to some 
interesting conclusions regarding the phenomenology of a 
compact star in scalar-tensor theory. 

The hrst thing that we should examine is the leading 
order behaviour of the orbital frequency of circular orbits 
and how this compares to the behaviour in GR. One can see 
that to leading order, the orbital frequency behaves as 


^ M- 


(54) 


This is the frequency one would measure for an object orbit¬ 
ing far from the central compact object, i.e. the “Keplerian” 
orbital frequency. So, if we were to measure the orbital mo¬ 
tion of matter orbiting far from the compact object in order 
to estimate its mass, the mass that we would infer assum¬ 
ing “Keplerian” motion would be M = M — Woao- It is 
the quantity M therefore that is the observable mass of the 
scalar-tensor theory and it is the mass of the compact object 
in the same sense as we would infer the mass in GR. 

Starting from that, we can attempt to compare the ex¬ 
pressions for AE, i}p/Q and 0,^/0, order by order in U and 


see if we can distinguish the orbital behaviour of particles 
between scalar-tensor gravity and GR. The hrst three terms 
of AE in GR read, 

+ (55) 

where U = with M being the mass as infered 

from Keplerian orbits in GR. The corresponding terms in 
scalar-tensor theory, as we have seen in Section 4, are 

-l_ (774 20Ji 5 

^3 V 9A72 9M 2 ) 9PP 

+0 (17®) , (56) 

where U = for the reasons discussed above. The 

comparison between the two expressions shows that, even 
when we can’t tell observationally if the mass that we are 
measuring corresponds to the GR M or the scalar-tensor 
M, we can still distinguish between GR and scalar-tensor 
theory. If we adopt the general form for the expansion of 
AE with respect to U to be 


= (57) 

k=2 


then the ratio 

A 4 f -3/2 , GR 

> scalar-tensor 

could be used as a tool to constrain deviations from GR 
or measure the scalar charge of a compact object. Gurrent 
constrains from solar system experiments and pulsar timing 
require that qq ^ 10“®. However, even for ao = 0 there is 
a discrepancy between GR and scalar-tensor theory that is 
equal to 2/3 o4Kq/ 3M^. Interestingly enough this deviation 
from GR manifests itself at an order where the rotation has 
no contribution yet. Therefore it is present for both rotat¬ 
ing and non-rotating stars and depends only on the scalar 
monopole and the parameter fio of the theory. To give a first 
estimate of the magnitude of the effect, we present in Ta¬ 
ble [T] the relative difference of the ratio A 4 /A 2 between GR 
and scalar-tensor theory, A(A 4 /A 2 ), dehned as the differ¬ 
ence of the GR value minus the scalar-tensor value divided 
by the GR value. The scalarized neutron star models used 
for the table are th e non-rotating models in the work by 
iDoneva et al.l |20l3). They are generated using the APR 
equation of state (EOS) for /3o = —4.5 and ifiao = 0. We 
can see that A(A 4 /A 2 ) is in the range of 1 % — 10 %, and we 
expect it to be even larger for rotating neutron stars. 

The same analysis can be performed for the expansion 
of Qp/fl. The first terms of the expansion in GR read 

+ (59) 

while the corresponding terms in scalar-tensor theory are, 
flp _ Wo [PoWo - 8aoM)\ ^ 4Ji 3 p/p4x 

2m5 

(60) 

One can see that the periastron precession at leading order 
behaves as flp/fl ~ 3U^ in GR while in scalar-tensor this 
is modihed. If we assume the expansion for the periastron 
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Table 1. Values for the mass M, the scalar charge Wq and the 
resulting deviations from GR, for non-rotating neutron star mod- 
els constructed with the APR EOS. The values were taken from 
iDoneva et al. Ilioii) and the models assumed /Sq = —4.5, ap = 0 
and a value for the scalar field at infinity, 0oo = 0. 


M/Mq 

M 

(km) 

Wo 

(km) 

A(A4/A2) 
(per cent) 

AB 2 

(per cent) 

1.68 

2.52 

-0.3596 

-4.07257 

-0.509071 

1.77 

2.655 

-0.4868 

-6.7236 

-0.84045 

1.86 

2.79 

-0.5448 

-7.62598 

-0.953248 

1.94 

2.91 

-0.5415 

-6.92534 

-0.865667 

2 

3. 

-0.4608 

-4.71859 

-0.589824 

2.05 

3.075 

-0.2414 

-1.23258 

-0.154072 


precession, 


np/n = Y,BkU\ (61) 

k=2 


then the relative difference between GR and scalar-tensor 
theory at leading order will be, 




Wo (/JqWo — SaoM'j 

6M2 


(62) 


Again if we impose the observational constrains that rednce 
Qo to be almost zero, we can see that the discrepancy be¬ 
tween GR and scalar-tensor theory remains (although it is 
smaller than in the case of the ratio (A4/A2) by a factor of 
8). Using the same numerical examples as before, we esti¬ 
mate in Table [ 1 ] the magnitude of the relative difference in 
B 2 between GR and scalar-tensor theory. 

The case of the nodal precession is slightly different, be¬ 
cause of the strict rotation dependence of the effect. The dif¬ 
ferences between GR and scalar-tensor theory comes about 
at higher order with respect to the multipole moments. In 
GR the nodal precession expansion has the form. 


n 


— U^ + 
M2 ^ 


3 M 2 
2M3 


U'^ + O ([/®) , 


(63) 


where we should note that the next term after the term of 
order in the expansion, is of order U®. In contrast, the 
nodal precession frequency in scalar-tensor theory has the 
expansion 


flz _2Ji 3 3(M2 — aoIU 2 ) 

IT ^ 2M^ 

_2.hWo{PoWyaoM)^,^^^^,^^ (64) 

where apart from the differences with respect to the U'^ 
term, there is also a term of order [7®. If we impose the 
observetional constrains for ao in this case, we can see that 
the terms of order and U* are the same between GR and 
scalar-tensor theory, but the U® term is non-zero in scalar- 
tensor. Unfortunately, this discrepancy probably comes at 
a high enough order to make measuring it in observations 
quite challenging. 


6 CONCLUSIONS AND OUTLOOK 

We have used the recently developed multipole moment for¬ 
malism in scalar tensor theory in order to obtain expressions 


for various observables that characterise geodesics in terms 
of the moments. Ryan had derived similar expressions in 
GR with the intention to use them for extreme mass ratio 
inspirals (EMRIs) as a tool to measure the multipole mo¬ 
ments of the spacetime around a supermassive BHs from 
gravitational waves observations. It is questionable whether 
the expressions we derived here will be useful in this par¬ 
ticular scenario i n scalar tensor t h eory. This i s because no ¬ 
hair theorems by|H awkind ( 197^) ; iBekenstei^ 1995l l ; 

iMavo fc BekensteinI 1 1996l i ; Sotiriou fc FaraoidI ( 2012 ), sug¬ 
gest that BHs will not carry a scalar charge. Realistic astro- 
physical BHs might circumvent no-hair theorems by violat- 
ing one or mor e of their assumptions, see a recent review by 
ISotirioul (l2015l l for a discussion. For instance, the presence of 
matter around the BH can introduce a scalar charge (blac k 
hole scalarizat i on) as it was shown bv I Cardoso et all ll2013ll : 
ICardoso et al.l (l2013l ) . Whether such a charge would be large 
enough to lead to some measurable effect is not clear at this 
stage. 

In principle one could consider using the relations be¬ 
tween the observables and the multipole moments in EM¬ 
RIs as a direct test of the BH scalarization hypothesis. 
Something similar could be done with supermassive BHs 
or stelar mass BH s usin g QPOs. After the first discovery 
bv iGierlihski et ^ (l2008l l of QPOs in an active galaxy, sev¬ 
eral more source s have been verified [see for example the 
review by[G uptal ll2014l ')1. These QPO sources could be used 
to test whether supermassive BHs are scalarized by imple¬ 
menting the expressions for th e frequencies presente d here . 
Additionally to these sources, Ijohannsen fc Psaltid (1201 il l 
have argued that Sgr A* could be another promising source 
of QPOs, measured using very-long baseline interferome¬ 
try, that could be used for this sort of tests. Apart from 
supermassive BHs, QPOs are commonly o bserved in X - 
ray binaries t hat h ost stellar mass BHs (see iLam 3 1 I 2003 I 1 : 
Ivan der KlisI ll2006l l'l. Therefore, with improved future ob¬ 
servations it might be possible to test the BH scalarization 
scenario. 

Since it is still questionable whether BHs can actually 
carry any significant scalar charge in scalar-tensor theory, 
the primary class of systems in which one could apply the 
derived relations are those when the central object is a NS. 
Therefore NSs in LMXBs are probably the most promis¬ 
ing sources for implementing the expressions developed here. 
These are usually members of a binary system, which pro¬ 
vides the opportunity for an independent relatively accurate 
measurement of the Keplerian mass. This could make even 
more successful the implementation of the extended Ryan 
expressions, since one would have less fitting parameters (the 
mass would be independently known) if one were to attempt 
a fit of t he QPOs frequ encies along the lines that was pro¬ 
posed bv IPapp^ (l2012li (additionally one would need to fit 
lower order coefficients than those discussed in the latter 
work). 

Furthermore one could attempt to associate the accre¬ 
tion spectrum and the disc temperature distribution, to the 
expression for AE in hope of probing the deviation between 
scalar-tensor theory and GR using some sort of disc tomog¬ 
raphy or the iron line rever beration technic (see for example 
work bv I Jiang et al.l (l2014l i'). 

It should be noted that since we considered here only 
observables associated with geodesics of the metric, our ap- 
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proach does not actually account for the effect of a scalar- 
scalar interaction between the central objects and the orbit¬ 
ing matter. We have essentially worked under the assump¬ 
tion that the orbiting matter is not scalarized, otherwise it 
would not follow geodesics of the metric in the first place. 

Further analysis is needed to explore the possibilities 
for using the expressions we have derived in order to extract 
constrains for scalar-tensor theories as a deviation from GR. 
This would require detailed modelling of the different classes 
of sources that we have discussed here and an artful consid¬ 
eration of current and future observational uncertainties. 
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APPENDIX A: FULL EXPRESSIONS FOR THE 
ISCO 

Here we present the complete expansion for the ISCO radii 
and the frequencies discussed in Sec. ESi 


Risco —6A/ ( 1 + 


0.452733aoWo 0.544331Ji 


A1 Scalar moments scaling as, W 2 - 
W 4 ~ 


and 


Pisco —' 


2V6m(^ 
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